UNCLASSIFIED 


Defense  Technical  Information  Center 
Compilation  Part  Notice 

ADPO 11590 

TITLE:  The  Application  of  Orthonormal  Electromagnetic  Beams  to 
Characterizing  Complex  Media 

DISTRIBUTION:  Approved  for  public  release,  distribution  unlimited 


This  paper  is  part  of  the  following  report: 

TITLE:  International  Conference  on  Electromagnetics  of  Complex  Media 
[8th],  Held  in  Lisbon,  Portugal  on  27-29  September  2000.  Bianisotropics 
2000 

To  order  the  complete  compilation  report,  use:  ADA398724 

The  component  part  is  provided  here  to  allow  users  access  to  individually  authored  sections 
of  proceedings,  annals,  symposia,  etc.  However,  the  component  should  be  considered  within 
the  context  of  the  overall  compilation  report  and  not  as  a stand-alone  technical  report. 

The  following  component  part  numbers  comprise  the  compilation  report: 

ADP01 1588  thru  ADP01 1680 


UNCLASSIFIED 


11 


The  Application  of  Orthonormal  Electromagnetic 
Beams  to  Characterizing  Complex  Media 

G.  N.  Borzdov 


Department  of  Theoretical  Physics,  Belarusian  State  University 

Fr.  Skaryny  avenue  4,  220050  Minsk,  Belarus 

Fax:  + 375-17-226  0530;  email:  borzdov@phys.bsu.unibel.by 

Abstract 

We  present  an  approach  to  characterizing  complex  media,  based  on  the  use  of  new  families 
of  orthonormal  electromagnetic  beams.  Each  familly  consists  of  orthonormal  exact  solutions 
of  Maxwell’s  equations,  which  differ  fundamentally  from  the  well-known  approximate  solu- 
tions — the  Hermite-Gaussian  and  Laguerre-Gaussian  beams.  A promising  type  of  such 
orthonormal  beams — beams  defined  by  the  spherical  harmonics — is  discussed.  The  proposed 
approach  makes  it  possible  to  use  beams  focused  into  a small  spot  on  the  sample  surface. 


1.  Introduction 

The  free-space  techniques  for  characterizing  complex  media  are  based  on  the  use  of  the  plane- 
wave  approximation  of  the  incident  beam.  Computer  modelling  [1,  2,  3]  of  the  free-space  tech- 
niques [1,  2,  3,  4],  based  on  the  covariant  impedance  methods,  has  shown  that  these  techniques 
make  it  possible  to  extract  all  material  parameters  of  an  anisotropic,  chiral,  or  general  bian- 
isotropic  medium,  provided  that  the  reflection  and  transmission  coefficients  of  planar  samples 
under  normal  and  oblique  incidence  of  plane  harmonic  waves  are  measured  with  sufficient  ac- 
curacy. However,  this  requires  a rather  complicated  measurement  setup,  and  in  many  cases  the 
plane-wave  approximation  of  beams  in  use  proves  to  be  inadequate,  especially  for  thick  samples. 

The  technique  presented  in  Refs.  [5,  6]  makes  it  possible  to  compose  a set  of  orthonormal 
beams  in  a complex  medium  or  free  space,  normalized  to  the  energy  flux  through  a given  plane. 
They  can  be  used  to  generalize  the  free-space  techniques  [1,  2,  3,  4]  for  characterizing  complex 
media,  developed  for  the  case  of  plane  incident  waves,  to  the  case  of  incident  beams.  A promising 
type  of  such  orthonormal  beams — beams  defined  by  the  spherical  harmonics — is  introduced  in 
Ref.  [7].  In  this  paper,  we  discuss  the  properties  and  applications  of  these  beams  in  more  detail. 


2.  Beams  Defined  by  the  Spherical  Harmonics 

In  this  paper,  we  consider  electromagnetic  fields  in  free  space  of  the  form  [7] 

Wj(r,  t ) = e~iu}t  dip  Je  eir-k^Y/(6,  <p)v{0,  <p)W(9,  ip)  sin  OdO. 

They  are  defined  by  the  spherical  harmonics 

Yr{e,<p)=NlmP\m\ cos0)e™* 
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and  Pfn{ cos 6)  and  ji{kr)  are  the  spherical  Legendre  and  Bessel  functions.  The  spherical  har- 
monics Y™(6,p)  satisfy  the  relations 

{YrW)  = r dp  r y,m* (0, p)YlV1' (0,  p)  sin OdB 
Jo  Jo 

= SwSmm'.  (4) 

Hence,  for  the  beams  under  consideration  (see  also  Ref.  [8]),  Bu  is  a unit  sphere  (Bu  = S2),  the 
beam  manifold  B C Bu  is  its  zone  with  0 € [0i,  02]  and  P G [0, 27r];  and  dB  = sin OdOdp. 

To  compose  electromagnetic  beams  in  free  space,  it  is  convinient  to  set  W = col(E,  B)  = 
col(E,  H).  For  a time-harmonic  field,  the  component  S3  = 63  • S of  the  time  average  Poynting 
vector  S can  be  written  as 

S3  = • (E*  x H + E x H*)  = W+QW,  (5) 

lD7 r 

where 

0__£_(  0 ei  g>  e2  - e2  8)  ei 

167r  y e2  8 ei  — ei  <g>  e2  0 

and  (g>  is  the  tensor  product.  Therefore,  for  the  electromagnetic  beams  Wj,  the  condition 
(W||<5|Wj)  = Nq  is  in  fact  the  normalization  to  the  beam  energy  flux  Nq  through  the  plane 
Co  normal  to  q = 63: 

(Wj|Q|Wj)  = [ SzdaQ  = Nq.  (7) 

J CTQ 

Each  family  of  the  fields  under  consideration  is  described  by  functions  which  have  integral 
expansions  in  plane  waves  with  wave  normals  lying  in  the  same  given  solid  angle  fi.  In  particular, 
one  can  set  the  angular  spectrum  of  plane  waves  by 

k = k(0,  ( p)  = k[0'(0,  ip),  p'{6,  p)),  (8) 

where 

k = k /k  — sin0'(ei  cosc£>'  + e2  sine/?')  + e3cos0'.  (9) 

In  this  paper,  we  restrict  our  consideration  to  beams  with 

o'  = KO0,  p'  = p,  (10) 


where  kq  is  some  real  parameter;  0 < no  < 1.  Correspondingly,  to  set  the  beam  base,  it  is 
convinient  to  use  the  radial,  the  meridional,  and  the  azimuthal  basis  vectors 


er(6',p)  = sin0'(ei  cos  p + e2sin<£)  + e3  cos0',  (11) 

e o'(0',p)  = cos  0'(ei  cosy?  + e2  sin  p)  — ezsinO1,  (12) 

e<p{p)  = — ei  sin p + e2  cos  p.  (13) 

Let  us  set  two  amplitude  functions  by 

W<^>S(b)  = (*)•  <14' 

w<0'^-(b  )=(-;,)■  (i5) 


Since  the  beams  with  the  amplitude  function  W [Eq.  (14)]  are  composed  from  plane  waves  with 
the  meridional  orientation  of  E and  the  azimuthal  orientation  of  B,  they  will  be  referred  to  as 
Em  beams  or  Ba  beams.  Similarly,  the  amplitude  function  W [Eq.  (15)]  results  in  E\  beams 
or  Bm  beams.  The  field  vectors  of  Em  and  Ea  beams  are  related  by  the  duality  transformation 
E -t  B,  B — ^ — E. 
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3.  Orthonormal  Beams 


Let  us  first  consider  orthonormal  beams  with  the  angular  spectrum  ft  = 2tt,  i.e.,  the  superpo- 
sitions of  eigenwaves  propagating  into  a given  halfspace.  To  this  end,  let  us  set  6i  — 0 and 
02  = 7r/2  in  Eq.  (1),  and  k0  = 1 in  Eq.  (10).  In  this  case,  the  amplitude  functions  W(6,ip) 
for  Em  and  Ea  beams  are  given  by  Eqs.  (14)  and  (15)  with  O'  — 6 , and  the  orthonormalizing 
function  v = v(0,  ip)  reduces  to  a constant  [5].  The  beam  manifold  B = Sjj  is  the  northern 
hemisphere  Sf?  of  the  unit  sphere  S2.  This  results  in  two  different  sets  of  orthonormal  beams 
defined  by  the  spherical  harmonics  Yf  with  even  and  odd  j,  respectively.  However,  it  is  of 
value  to  have  a complete  system  of  orthonormal  beams  Wj  defined  by  the  whole  set  of  spherical 
harmonics  Y J,  for  which  (Wj|Q|W^)  = 0,  if  at  least  one  of  the  three  conditions  is  met:  j'^j, 
s'^s,  or  the  beams  have  the  alternative  polarization  states  (Em  and  Ea  beams). 

To  this  end,  let  us  set  the  beam  base  of  time-harmonic  electromagnetic  Em  beams  Wj  [Eq. 
(1)]  by  Eqs.  (9)  and  (10)-(15)  with  0\  = 0,  62  = 7r,  and  kq  < 1/2.  In  this  case,  the  beam 
manifold  is  the  unit  sphere*  (£?  = S2),  the  angular  spectrum  ft  < 27T,  and  the  orthonormalizing 
function  has  the  form  


__  2 I 2ttkqNq  sin(Ko^) 
A V c sin  6 


(16) 


These  Em  beams  also  can  be  expanded  into  a series  as  described  in  Ref.  [5].  As  before,  Em  and 
Ea  beams  are  related  by  the  duality  transformation. 

The  smaller  is  kq?  the  smaller  is  the  angular  spectrum  fi,  i.e.,  the  more  collimated  is  a beam. 
Conversely,  if  ko  = 1/2,  i.e.,  ft  = 2i r,  the  beam  becomes  highly  localized  and  has  an  energy 
distribution  in  the  core  region  similar  to  the  beams  presented  in  Refs.  [5,  6,  7].  When  0 and 
k = 1/2,  or  k ^1/2,  these  beams  resemble  electromagnetic  tornadoes  with  spiral  energy  fluxes 
and  pronounced  core  regions. 

The  general  time-harmonic  beam  with  two-dimensional  beam  manifold  B can  be  written  as 


W(r,i)  =e~iwt 


L 


eirk{b)v(b)u(b)W(b)dB , 


(17) 


where  u : B — > Cl  is  a complex  scalar  function  on  B . Let  (un)  be  an  orthonormal  base  of 
complex  functions  on  B.  Then,  the  function  u can  be  expanded  into  a series  as 


u(b)  ='^2  cnun(b), 

n 


(18) 


where  Cn  = (um\u).  By  applying  the  approach  described  in  Refs.  [5,  8],  we  obtain  an  expansion 
of  W (17)  into  a series  of  orthonormal  beams  Wn  as 

W = £c„Wn.  (19) 

n 


It  is  essential  that  the  coefficients  Cn  can  be  extracted  from  the  beam  W as  follows: 

= jL<W„|Q|W).  (20) 

What  is  even  more  important  they  are  measurable  values  provided  that  there  exists  a source 
of  orthonormal  beams  W„.  As  it  is  shown  in  Refs.  [5,  7],  X = (W|Q|W)  is  the  energy  flux 
through  the  plane  cto  in  the  case  of  time-harmonic  beams  with  two-dimensional  manifold  B. 
Each  of  the  complex  coefficients  Cn  of  the  beam  W (19)  can  be  calculated  from  the  results  of 
three  measurements  [5,  7]. 
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4.  Conclusion 

The  presented  approach  to  characterizing  complex  media  is  based  on  the  use  of  new  families 
of  orthonormal  electromagnetic  beams.  Each  familly  consists  of  orthonormal  exact  solutions 
of  Maxwell’s  equations,  which  differ  fundamentally  from  the  well-known  approximate  solutions 
— the  Hermite-Gaussian  and  Laguerre-Gaussian  beams.  By  using  these  solutions,  the  results 
obtained  in  Refs.  [1,  2,  3,  4]  for  the  case  of  plane  incident  waves  are  generalized  to  the  case 
of  time-harmonic  beams  obliquely  incident  onto  a general  bianisotropic  slab.  To  this  end,  the 
fields  of  incident,  reflected,  and  transmitted  waves  are  expanded  into  series  of  orthonormal  vector 
functions.  The  obtained  solutions  make  it  possible  to  calculate  the  complex  scalar  coefficients 
of  these  series.  It  is  shown  that  these  coefficients  are  measurable  values,  and  the  corresponding 
measurement  scheme  is  suggested.  Assuming  that  they  are  given  or  measured,  it  is  possible 
to  reconstruct  the  reflection  and  transmission  coefficients  of  the  slab  for  partial  incident  plane 
waves  and  then,  using  the  techniques  presented  in  Refs.  [1,  2,  3,  4],  to  extract  the  whole  set 
of  material  parameters.  One  can  use  various  families  of  orthonormal  beams,  in  particular,  the 
family  of  beams  defined  by  the  spherical  harmonics.  Results  of  numerical  analysis  of  the  latter 
and  peculiarities  of  its  possible  application  to  characterizing  various  complex  media  will  be 
reported  orally. 

The  proposed  approach  makes  it  possible  to  use  beams  with  wide  angular  spectrum,  focused 
into  a small  spot  on  the  sample  surface.  Usage  of  well  focused  beams  eliminates  the  need  to 
work  in  an  anechoic  enviroment. 
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